Phonon spectral function for an interacting electron-phonon system 
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Using exact diagonalzation techniques, we study a model of interacting electrons and phonons. 
The spectral width of the phonons is found to be reduced as the Coulomb interaction U is increased. 
For a system with two modes per site, we find a transfer of coupling strength from the upper to the 
lower mode. This transfer is reduced as U is increased. These results give a qualitative explanation of 
differences between Raman and photoemission estimates of the electron-phonon coupling constants 
for AaCeo (A= K, Rb). 



In a metallic system a phonon can decay into electron- 
hole pair excitations. This decay contributes to the 
width of the phonon. It was pointed out by Allen that 
this additional broadening, can be used to estimate the 
electron-phonon coupling.EI The width can be measured 
in neutron scattering or, for the orien±ationally disor- 
dered fuUerenes, in Raman scattering.B Normally, the 
electron-phonon couuling is deduced by assuming non- 
interacting electronsE The method is, however, often 
applied to systems with strong correlation due to the. 
Coulomb interaction, such as the alkali-doped fuUerenes .tl 
In the alkali-doped fuUerenes the electron-phonon inter- 
action plays an important role, and accurate estimates 
of the coupling strength are essential. Almost all exper- 
imental estimates for these systems are based on Allen's 
formula, and the accuracy of this formula is therefore 
crucial. 

In strongly correlated systems the hopping is reduced 
and the excitation of electron-hole pairs may be more dif- 
ficult. For instance if the correlation is so strong that the 
system has a metal-insulator transition, the decay into 
electron-hole pair excitations is completely suppressed. 
One aim of this paper is therefore to study how the esti- 
mate of the electron-phonon coupling is influenced if the 
electron-electron interaction is taken into account. 

In metals a phonon can decay into a (virtual) electron- 
hole pair excitations which can then decay into a dif- 
ferent phonon. In this way there is a coupling between 
different phonon modes of the same symmetry, leading 
to new modes which are linear combinations of the old 
ones. These new modes can have quite different coupling 
strengths than the old modes. A second aim of this pa- 
per is to study how the coupling strength is transferred 
between the modes due to the coupling via electron-hole 
pair excitations. 

The electron-phonon coupling has been studied ex- 
tensively for the alkali-doped fuUerenes. In particulat, 
there have been a study based on neutron scattfijingjCl 
and several studies based on Raman scatteringJIij The 
high resolution studies of Winter and KuzmanyH show a 
very strong coupling to a few of the low- lying modes, 
but almost no coupling to the high-lying modes. An 
alternative approach is based on photoemission from 
free negatively charged C^q molecules. By studying the 
weight of vibration satellites, it is possible to deduce the 



electron-phonon coupling.Q This results in rather differ- 
ent electron-phonon coupling constants. Although the 
main coupling was to the low-lying modes, there was also 
a substantial coupling to the two highest Hg modes. The 
total coupling strength was also larger than deduced from 
Raman scattering. 

In this paper we study a simple model with electron- 
phonon and electron-electron interactions. We consider 
a finite cluster with a nondegenerate electronic level and 
a nondegenerate phonon on each site. This model is 
solved by using exact diagonalization. We find that the 
Coulomb interaction reduces the phonon width, and that 
the use of Allen's formulatl therefore leads to an underes- 
timate of the electron-phonon coupling constants in Ra- 
man scattering experiments. Furthermore we find that 
due to the indirect interaction of different phonon modes 
via electron-hole pair excitations in metallic systems, 
there is a transfer of coupling strength to the low-lying 
modes which is not present for a free molecule. Since the 
Raman measurements of the electron-phonon coupling 
are performed for a solid, but the photoemission estimate 
is for a free molecule, the weight transfer is present in the 
Raman but not in the photoemission estimate. These 
observations are consistent with differences between the 
coupling constants deduced from Raman scattering and 
photoemission. 



We consider a model with Nn 
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nondegenerate 



phonons per site and with electrons without orbital de- 
generacy. The Hamiltonian is 
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where i labels the A^sitc sites, Cia- and bi^ annihilate an 
electron with spin a and a phonon with the label i^, re- 
spectively, on site i and rii^ = H^^ia is an occupation 
number operator. The energy of the phonon v is w^ and 
its coupling to the electrons is is 5^. The correspond- 
ing dimensionless electron-phonon coupling is given by 
A^ = 2glN{0)/ujph, where N{0) is the density of states 
per spin. The energy of the electronic level is cq. Two 
electrons on the same site have a Coulomb repulsion U. 
The hopping between the sites is described by matrix 



elements tij. A Hamiltonian like (P with tij = t 



for 



the nearest neighbor hopping has a high symmetry and 
a correspondingly large degeneracy. Since we use exact 
diagonalization to solve the model, we have to limit the 
number of sites to a small number (4-6). The result- 
ing one-particle states are then very sparse in energy. 
Therefore we lower the symmetry by choosing each tij 
randomly within some interval, which leads to a denser 
energy spectrum. The model is solved and the result is 
then averaged over different sets of {Uj}. The strength 
of the hopping is measured by the one-particle width 
W of the electronic band (for g — and [/ = 0). In 
this model we for simplicity consider nondegenerate (Ag) 
phonons and electrons. In, for instance, A3C60 (A= K, 
Rb) the phonons are five-fold degenerate Hg phonons and 
the electron states have a three-fold orbital degeneracy. 
The new physics which may be introduced by these de- 
generacies, e.g., the Jahn- Teller effect, is not considered 
here. 

We consider a half- filled system, i.e., A'sitc electrons. 
With A'sito — 6 there are then 400 different electronic con- 
figurations. To obtain a finite size Hilbert space, we limit 
the maximum number of phonons per modes to A^phon- 
The number of phonon states is then (A^phon + 1) ='"= 
for the case of one mode per site. For instance, with 
A^gito — 6 and A^phon — 3, the total Hilbert space has the 
dimension 1.6384 • 10^. Such a problem can be solved 
using exact diagonalization, i.e., the ground-state is ex- 
pressed as a linear combination of all possible basis states 
in the Hilbert space. The lowest eigenfunction of the 
corresponding Hamiltonian matrix is then found using 
the Lanczos method. We further calculate the phonon 
Green's function 



figure illustrates that one underestimates the electron- 
phonon coupling if Allen's formula is used to extract the 
coupling for a system with a finite U. For systems like 
A3 Ceo (A= K, Rb), where the Coulomb interaction is be- 
lieved to play an important role, the width of the phonons 
may then be substantially reduced. The use of Allen's 
formula would then correspondingly underestimate the 
electron-phonon coupling. 




FIG. 1. The phonon spectral function for different values 
of the Coulomb interaction U. The insert shows the width 7 
of the spectral function as a function of U. The figure illus- 
trates the narrowing of the spectral function as U is increased. 
The parameters are uji — 1, X — 0.073 and W = 2.5 and the 
system has six sites. 
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where |0) is the ground state, 0i,y(i) — hiy{t) + h\^{t) is 
the phonon field operator in the interaction representa- 
tion and T is the time-ordering operator. Calculation 
of the Fourier transform gives D^{ijj). We then define a 
spectral function as 
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and study the average pp,i(w) = X],yi ^M«''(^)/^si*e- Due 
to the finite size of the system, the spectrum is discrete. 
We therefore introduce a Lorentzian broadening with the 
FWHM (full width at half maximum) 0.01 eV. 

Fig. shows the phonon spectral function A{ijj) for 
different values of U . Due to the small size of the sys- 
tem, the width of the spectrum should not necessarily be 
expected to agree with Allen's formula even for U — Q. 
Nevertheless, the result of AUen's formula 7AUcn = 0.19, 
is comparable to the width found for f/ = 0. The figure 
illustrates how the spectral function becomes narrower 
with increasing U . This is further illustrated by the in- 
set, which shows the width of the spectrum, calculated as 
the mean square deviation of the spectral function. The 



We next discuss the case when there are two phonon 
modes per site, which have the unperturbed energies u)i 
and ijJ2- First we calculate the lowest order phonon self- 
energy. This involves evaluating a "bubble" diagram. 
The self-energy can be written as 
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where /(w) depends on the precise band structure. We 
consider contributions to the self-energy which are both 
diagonal and non-diagonal in the index ly. The non- 
diagonal contribution corresponds to a phonon v decay- 
ing into an electron-hole pair followed by this electron- 
hole pair decaying into a phonon v . The non-interacting 
phonon Green's function is 
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The interacting phonon Green's function is then given by 
D-\^)^{D\l.)\-^-JI{^)^ (6) 
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The modes of the coupled system are obtained by looking 
for zeros of the determinant of the matrix in Eq. (^ . For 



the lowest mode, the corresponding eigenvector consists 
of a bonding Hnear combination of the two unperturbed 
modes. As a result the coupling to the electrons is in- 
creased for these modes, due to constructive interference 
between the couplings for the two unperturbed modes. 
In the same way the coupling is reduced for the higher 
mode. For instance, we can look for the width of the 
lowest mode in the limit when lo^ ^ wi and when the 
electron-phonon coupling is weak. We then find that the 
width of the lowest mode is increased by a factor of 

(I + CA2), (7) 

and the width of the highest mode is reduced by a factor 

(8) 
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where c is somewhat larger than unity (c ^ 3) and de- 
pends on the shape of the band. 




FIG. 2. The phonon spectral function for a system with 
one phonon mode (upper part) and two phonon modes (lower 
part) per site for [7 = 0. The insert shows the width of the 
lower mode as a function of U in the cases of one or two 
modes per site. Comparison of the widths in the two cases, 
illustrates how the lower mode is broadened for small values 
of f/ due to the interaction with the upper mode. The param- 
eters are W = 3.7, oji = 0.5, uj2 — 1, gi — 0.3 and 32 = 0.4 
and the system has four sites. AH spectra have been given a 
Lorentzian broadening with the FWHM=0.1 

The result in Eq. (R) is based on the lowest order 
phonon self-energy and it neglects the Coulomb repul- 
sion completely. We therefore study the same problem 
using exact diagonalization. Fig. ^ compares results for 
systems with one or two modes per site. The discrete 
spectra have been broadened by a Lorentzian with the 
FWHM=0.01. The main figure shows results for [/ = 0, 
and it illustrates how the lower mode is broadened when 
the higher mode is switched on. For the parameters in 
Fig. I Ai = 0.043 and A2 = 0.085, and the additional 



broadening of the lowest mode is of the order of magni- 
tude predicted by Eq. (|7|). The insert shows the width 
of the lower mode as a function of U. These results were 
obtained by fitting Lorentzians to the broadened spectra. 
The width for very large values of U is due to the broad- 
ening of the discrete spectrum that we have introduced. 
As U is increased, the width of the mode is reduced, as 
discussed above. The figure further illustrates that the 
transfer of coupling strength is reduced as U is increased. 
This is expected, since the effects of hopping, and thereby 
the indirect coupling, is reduced as U is increased. 

It would be interesting to repeat these calculations for 
systems with degenerate phonons, e.g., to include the 
Jahn- Teller effect. This leads, however, to systems which 
are so large that they cannot easily be treated using exact 
diagonalization. Within a Hartree calculation we find a 
similar transfer of coupling strength to the lower modes 
also for Jahn- Teller phonons and electrons with orbital 
degeneracy. The transfer is, however, reduced by the 
nonspherical parts of the Coulomb interaction, i.e., by 
the difference between the interaction for equal orbitals 
and unequal orbitals. This effect may also play a role 
when we go beyond the Hartree approximation. 

Finally, we observe that in theoretical approaches 
which do not explicitly include the transfer of coupling 
strength between the modes, it is appropriate to include 
this transfer by using the corresponding coupling con- 
stants. On the other hand, in a treatment where this 
transfer is explicitly included, the transfer should not be 
contained in the coupling constants used in the model. 

To summarize, we have calculated the phonon spec- 
tral functions for systems with interacting electrons and 
phonons. We find that the Coulomb interaction between 
the electrons reduces the width of the phonons caused 
by the phonon decay into electron-hole pairs. As a re- 
sult, estimates of the electron-phonon coupling based on 
the phonon width underestimate this coupling unless the 
Coulomb interaction is taken into account. This is con- 
sistent with the observations that weaker couplings have 
been deduced from Raman measurements than from pho- 
toemission (PES) experiments. Furthermore, we find 
that there is a transfer of coupling strength from the 
higher modes to the lower modes due to an indirect 
interaction via electron-hole pairs. This may, at least 
partly, explain the difference in the distribution of cou- 
pling strength between Raman and PES estimates, al- 
though it can probably not fully explain the weak cou- 
pling to the two highest phonons seen in Raman spec- 
troscopy. In this work we have treated nondegenerate 
phonons. It would be interesting to extend the work to 
degenerate, Jahn- Teller phonons, since these are the im- 
portant phonons in the alkali-doped Fullerenes. 
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